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0^ 

Abstract. To each representation of the elliptic quantum group Er,r)[sl2) 
is associated a family of commuting transfer matrices. We give common 
eigenvectors by a version of the algebraic Bethe ansatz method. Special 
cases of this construction give eigenvectors for IRF models, for the eight- 
vertex model and for the two-body Ruijsenaars operator. The latter is 
, a g-deformation of Hermite's solution of the Lame equation. 

1. Introduction 

> , 

'nI" , The Bethe ansatz is a method to construct common eigenvectors of com- 

I muting famihes of operators (transfer matrices) occurring in two-dimensional 

ly-^ ■ models of statistical mechanics. Faddeev and Takhtadzhan |TF] reformu- 

\ lated the problem into a question of representation theory: commuting fami- 

' lies of transfer matrices are associated to representations of certain algebras 

with quadratic relations (now called quantum groups). Eigenvectors are 
constructed by properly acting with algebra elements on "highest weight 
■ vectors". In this form, the Bethe ansatz is called algebraic Bethe ansatz. 

Q^, Whereas this construction has been very successful in rational and trigono- 

metric integrable models, its extension to elliptic models has been problem- 
atic, although the Bethe ansatz, in the case of the eight-vertex model, is 



> 

X 



^ , known since Baxter's work ||B|]. For elliptic models associated to s/at, an al- 

d ' gebra with quadratic relations has been introduced by Sklyanin , but the 

notion of highest weight vector is not defined for its representations, which 
makes a direct application of the algebraic Bethe ansatz impossible. 

Recently, a definition of elliptic quantum groups Er^nid ) associated to any 
simple classical Lie algebra q was given |^] . It is related to a g-deformation of 
the Knizhnik-Zamolodchikov-Bernard equation on tori. The representation 



theory of -Er.r? (■5/2) was described in | FV3 | 



In this paper we describe the algebraic Bethe ansatz for Er^n{sl2)- The 
construction is very close to the construction done in the trigonometric case 



in |TF|. The main difference is that transfer matrices act on spaces of 
vector-valued functions rather than on finite dimensional vector spaces. 

The basic results in this paper are Theorems ^ and ||. We present them 
in Sect. ^ after a summary of the representation theory of Er^-qish) in Sect. 
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^. In Sect. 1^ we formulate our result in the discrete case, which gives a 
construction of eigenfunctions for interaction-round-a-face models (Sect. ||). 

In Sect. ^ we introduce a version of Baxter's vertex-IRF transformation. 
With the help of this, we show how to obtain eigenvectors of the eight- 
vertex model transfer matrix from our eigenvectors. From this point of 
view, our construction is very reminiscent to Baxter's transformation of the 
eight-vertex model into an inhomogeneous six- vertex model [Q. Also, this 
vertex-IRF transformation shows a very close relation of E^.r^ with Sklyanin's 
algebra in the s/2 case. Note however that Sklyanin's algebra can only be 
defined for sl]\f, whereas the elliptic quantum group exists (at least) for all 
classical simple Lie algebras. 

Another class of problems in which the Bethe ansatz has been applied 
successfully is the class of Calogero-Moser-Sutherland quantum many body 
problems on the line with elliptic potentials, see [ FV1 |, [FV2]. In the case 
of two bodies, the Bethe ansatz goes actually back to Hermite, who solved 
in this way the generalized Lame equation, cf. [WW]. These integrable 
Schrodinger operators admit a g-deformation due to Rujsenaars In 
Sect. we present a (7-deformation of Hermite's result, i.e., we give eigen- 
functions for the two-body Rujsenaars operator. These eigenfunctions are 
parametrized by a spectral curve, similarly to the differential case. It is a 
double covering of a hyperelliptic curve. This result follows from our result 
and the observation that the transfer matrix associated to an evaluation 
representation is equal to the Ruijsenaars operator up to a scalar factor de- 
pending on the spectral parameter. A similar observation, relating transfer 
matrices to Ruijsenaars operators, has been made recently by Hasegawa Q 
in the context of Sklyanin's algebra, indicating that our construction should 
be extendable to the A^-body case. 



2. Modules over Er,r){sh) and transfer matrices 

In this paper, we construct eigenvectors of the transfer matrix of the ellipic 
quantum group Er^-qish)-, 0, FV5 ], associated to certain highest weight 
modules. 

Let us recall the definitions: we fix two complex parameters r, r/, such 
that Im(r) > 0. The definition of Er^nish) is based on an i?-matrix R{z, A) 
which we now introduce. Let 



^^i(j+l)V+27ri(j + i){t+i) 



be Jacobi's theta function and 



g(A + 2r?)g(z) eiX + z)9{2rj) 

Let y be a two dimensional complex vector space with basis e[l],e[— 1], 
and let Eije[k] = 5jke[i], h = En — Then, for z,X G C, R{z,X) G 
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End(V V) is the matrix 

R{z,X) = E^ii + _i _i + a(z, A)£;ii _i 

+ a{z, -A)^_i _i <g) Ell + fi{z, \)Ei^^i ® E^i^i + (3{z, -\)E^i^i ® E^i _i. 

This i?-matrix obeys the dynamical (or modified) quantum Yang-Baxter 
equation 

_ A- 2r//i(3) ) i2(i3) (^^ A) i?(23) A- 2j?/i(i) ) = 

in End(y(E)V(E'y), z, it;, A S C. The meaning of this notation is the fohowing: 
i?(^^)(A — 2rih^^y)vi ^V2 0vs is defined as 

{R{z, A - 2r]fi3)vi (g) V2) (8) V3, 

if hv3 = /is^s- The other terms are defined similarly: in general, let 
Vi, . . . ,Vn be modules over the one dimensional Lie algebra f) = Ch with 
one generator h, such that, for all i, Vi is the direct sum of finite di- 
mensional eigenspaces Vi[ij] of h, labeled by the eigenvalue /x. We call 
such modules diagonalizable P) -modules. If X G End(Fi) we denote by 
G End(Vi • • • K) the operator • • • Id (g) X Id (g) • • • acting 
non-trivially on the ith factor, and if X = Yl-^k 'S' £ End(Vi Vj) 
we set ,fj,n) is a function with values in 

End(yi • • • K), then X{h(^\. . . , = . . . , ;U„)i; if /i^t; = ^liV, 

for all i = 1, . . . , n. 

Now, by definition, a module over E-r^rjish) is a diagonalizable f) -module 
W = e^xecW^N, together with an L-operator L{z,X) € End|, ® W)) 
(a linear map commuting with h^^^ + h^"^^) depending meromorphically on 
z, A G C and obeying the relations 

{z-w,X- 27]h'^^^ ) (z. A) L(23) (z^;, A - 2r//i(i) ) = 

For example, W = V, L{w,X) = R{w — zo,X) is a module over Er^nish), 
called the fundamental representation, with evaluation point zq. In [[FV3f| 
more general examples of such modules were constructed: in particular, 
for any pair of complex numbers A, z we have an evaluation Verma module 
Va(z). Also, we have a notion of tensor products of modules over Er^nish)- 
The main examples considered in this paper will be tensor products (zi)® 
• • • (g V\^{zn) of evaluation Verma modules and some of their subquotients. 

For any module W over E-r^nish), we define the associated operator alge- 
bra, an algebra of operators on the space Fun(Ty) of meromorphic functions 
of A G C with values in W. It is generated by h, acting on the values, and 
operators a{z) , b{z) , c{z) , d{z) . Namely, let L{z) G End(y Fun(VF)) be the 
operator defined by {L{z){v /))(A) = L{z, X){v /(A — 2r]ii)) if hv = fiv. 
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View L{z) as a 2 by 2 matrix with entries in End(Pun(VF)): 

L{z){e[l]^f) = e[l](S)a{z)f + e[-l](S)c{z)f, 
L(z)(e[-1]®/) = e[l]<S,biz)f + e[-l]<S)d{z)f. 
The relations obeyed by these operators are described in detail in [ FV3[ | (in 



[ FVS ] these operators are denoted by d{z),b{z) and so on). 

Theorem 1. [^V|l For any module W over -E'r,r?('5/2); the transfer matri- 
ces T{z) = a{z) + d{z) preserve the space H = Fun(VF)[0] of functions 
with values in the zero weight space W[0], and commute pairwise on H: 
T{z)T{w) = T{w)T{z) on H. 



Theorem 2. |FV3| Let W = Va^(2;i) (8) • • • Va^(z„) he a tensor product of 
evaluation modules, and /et A = Ai + • • • + A„. Then W[K\ = Cfo, where vq, 
viewed as a constant function in Pun(VF) obeys the following highest weight 
condition: for every z, c{z)vo = 0,a{z)vo = A{z, X)vo,d{z)vQ = D{z,X)vq, 
with highest weight functions 

fj^ ^(^ - Qj) 

where we set pj = zj + ?/(— Aj + 1), qj = Zj + r?(Aj + 1) 

A highest weight module W of highest weight {A,A,D), is a module 
with a highest weight vector vq G Fun{W) such that c{z)vo = 0, a{z)vo = 
A{z, X)vo, d{z)vo = D{z,\)vo, for all z,X, and so that Fun(M^) is spanned 
by the vectors of the form • • • b{tj)vo, as a vector space over the field of 



meromorphic functions of A. It is shown in |FV3| that if A, D are of the form 
(^) , for some Pk,Qkj then every irreducible highest weight module of weight 
(A, A, D) is isomorphic to a subquotient of (8>VAj(zi), where the parameters 
Zj,Aj are related to pk, Qk as in the theorem. If p/., qi. are generic under the 
condition that Y2iPi ~ Qi) — ^^/A, then all highest weight modules of weight 
{A,A,D) are isomorphic to '^V\.{zi) itself. 

3. Bethe ANSATZ 

In this section we fix a highest weight module W of weight (A, A, D) of 
the form with highest weight vector vq. We assume that A is an even 
integer 2m > 0, so that the zero- weight space W[Q] can be nontrivial. 



We follow the strategy of [rF|: we seek common eigenvectors of T{uj) 
in the form b{ti) ■ ■ ■ b{tm)v, where v € Fun(VF)[A]. The problem is to find 
conditions for ti, . . . ,tn,v so that we have an eigenvector. The question of 
completeness, i.e., whether "all" eigenvectors can be obtained in this way, 
will not be addressed here, except in the example of the g-analogue of the 
Lame equation discussed below. 

Any non-zero vector v £ Fun(l^)[A] is of the form v = g{X)vQ, for some 
meromorphic function g ^ 0- To find for which values oi ti, . . . ,tm and which 
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choice of v we get an eigenvector we commute a and d with the 6's using 
the relations of the quantum group. The relevant commutation relations are 
(see [|V|) 

a{w)b{t) = r{t - w,X)b{t)a{w) + s{t - w,X)b{w)a{t), 

d{w)b{t) = r{w -t,X-27]h)b{t)d{w) - s{t-w,X-2rih)b{w)d{t). 

The functions r, s are defined by the formulae 

_ 9{t - 2r,)e{X) _ e{t + X)e{2r]) 

^ ' ^ " 9{t)eiX - 2r?) ' ' ~ 9{t)9iX - 2r?) ' 

Note that these coefficients have poles for t € Z + rZ, so we have to be 
careful to avoid infinities in the calculations: we assume that w,ti, . . . ,tfn 
are all distinct modulo Z + rZ. Using repeatedly the commutation relations 
to bring a to the right of the 6's, we get 

a{'w)b{ti) ■ ■ ■ b{tm) = Aobiti) ■ ■ ■ b{tm)a{w) 

n 

+ Ajb{ti) ■ ■ ■ b{tj.i)b{w)b{tj+i) ■ ■ ■ b{tm)a{tj) 
i=i 

for some complex coefficients Aj = Aj{w,ti, . . . ,tm, X). The first term is 
called the "wanted" term, and the others are "unwanted" terms. Similarly, 
we have 

d{w)b{ti)---b{tm) = Dob{ti)---b{trn)d{w) 

n 

+ ^ D^bitx) ■ ■ ■ b{tj-i)b{w)b{tj+i) ■ ■ ■ b{tm)d{tj) 

for some complex coefhcients Dj = Dj{w, ti, . . . , tm, A). The coefficients Aq 
and Ai are easy to compute, since they are products of coefficients appearing 
in the commutation relations: 

m 

Ao = llr{tj-w,X + 2r,{j -1)), 

m 

Ai = s{ti-w,X)llr{tj-ti,X + 2r]{j -1)). 
i=2 

A similar calculation gives 

m 

Do = '[[r{w-tj,X-2r,{j -1)), 

m 

Di = -s{ti-w,X)Wr{ti-tj,X-27]{j -I)). 
i=2 

A direct calculation of the other coefficients is more complicated. However 
the answer is simple, thanks to the 
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Lemma 3. For any permutation a of m letters, and all j = 1, . . . ,m, 

^j(w^)*(T(l), • • • )i<T(m)) A) = A^(^j){w,ti, . . . ,tm,X), 
• • • )*cr(m)i = D^(^j-j{w,ti, . . . ,tm,X), 

The proof of this Lemma is deferred to Sect. |8|. 

The next step is to find conditions for cancellation of unwanted terms. 
We have, with v{X) = g{X)vo, 

T{w)b{ti)---b{t^)v = Cob{ti)---b{tm)v (3) 

n 

+ Cjb{ti) ■ ■ ■ b{tj-i)b{w)b{tj+i) ■ ■ ■ b{tm)v 

i=i 

for some coefficients Cj. The condition of cancellation is Cj = 0, j > 1. Let 
us first consider Ci. It is given by 

The condition Ci = is then equivalent to 

r{t2-ti,X + 2?7) • • • r{tm-ti,X + 2r/(m - l))g{X + 2r/(m - 1)) 
= r{ti -t2, A - 2?7) • • • r{ti -t^, X - 2r?(m - l))D{ti,X + 2i]m)g{X + 2r/(m + 1)) 

which can be written as 

1^ ejtj -ti- 27?) -A- ejti - Qk) 9{X + 2i]{m - l))e{X + 2r,m)g{X + 2r]{m - 1)) _ 

|i ^^^^ - *i + 2^?) feii - P'^^ ^(^)^(^ - + ^^^'^ + 

In particular, the left-hand side of this equation should be independent of 
A. This holds if g is taken in the form 

^ - 0(A-2r?i) 



,=1 ^(2^) 

Thus C\ vanishes if g has this form and the ti obey the equation 

n ejtj -ti- 2r?) -A- 6l(ti - Qk) ^ 

Using Lemma ^ we find the conditions for Cj to vanish for j = 1, . . . ,m. 
The result is: 

Theorem 4. Lei W be a highest weight module over Er^nish) of highest 
weight {A,A,D) with A,D of the form @ and A = 2m € 2Z>o. Let 
T{w) G End(if), w £ C, be the corresponding transfer matrices. Let 
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v{X) = e'^^Y[T=i^i^ ~ '^'nj)vo- Then, for any solution {ti, . . . ,tm) of the 



Bethe ansatz equations 



lle{t,-t, + 2v)l\9{t.-Pk)~ ' w 

such that, for all i < j, ti ^ tj mod Z + rZ, the vector b{ti) ■ ■ ■ b{tm)v is a 
common eigenvector of all transfer matrices T{w) with eigenvalues 



2vc TT -w-2r]) TT ^ih - w + 2rj) yr e{w-pk) 

U e{t,-w) + n e{t,-w) l}^eiw-qk) 



What is left to prove is the formula for the eigenvalue, which is given by 
Co in (^). By definition 

^ g(A + 27?(m-l)) , ^ g(A + 2r?(m + 1)) 
= ^° g(A + 2r,m) + g(A + 2r,m) ^ + 

By inserting the formulas for Aq, Dq, g, D, we see that the A dependence 
disappears for the same reason as before, and we are left with the formula 
for e given in the theorem. 

We conclude this section by giving an explicit formula for 6(ti) * * * b(t^)v. 
It is sufficient to consider the case where W = Vai{zi) • • • V\„{zn), 
for any other highest weight module considered in the previous Theorem 
is a quotient of the submodule of a module of this form generated by the 
product of highest weight vectors. Recall |FV3] that Va{z) is defined to 
be the infinite dimensional vector space with basis cq , ei , 62 , . . . , such that 
hcj = (A — 2j)ej, and with the action of the generators of Et-^j^[s12) given 
by explicit formulae. 

Theorem 5. With the notations as in the previous theorem, 

b{t,)...b{t^)v = (-ire^(A+2.™) E nn n 

h,...,I„ 1=1 i£h k=l+l ^ * '^^^ 

^ n n ~ ~ ^''^ 



n 

X 



nn 



0{ti - tj) 

e{\ + tj-qk + 2r]mk-2r]YJI=k+i{^l-'^^l)) 



Of 4. \ ^nii'^ 

The summation is over all partitions of {1, . . . ,m} into n disjoint subsets 
/i ,...,/„ . The cardinality of Ij is denoted by rrij . 

In particular, we see that our eigenvector ^(A) is an entire function of A 
obeying + 1) = {—l)'^e^ip{\). The proof of this Theorem is contained 
in Section S. 
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4. Discrete models 

The construction of the transfer matrix admits the following variation. 
The difference operators a{w), . . . ,d{w) shift the argument of functions by 
ib2r/. Therefore we may replace Fun(VF) by the space Fun^(Ty) of all func- 
tions from the set = {// + 2T]j,j G Z} to W. The operators are then 
well-defined on Fun^(Ty) if jj, is generic. Also, it follows from Theorem ^ 
that the restriction to of the Bethe ansatz eigenfunctions is well defined 
for all fi. We thus have: 

Corollary 6. Suppose ti, . . . ,tm is a solution to the Bethe ansatz equations 
(0). Then, for generic fj,, the restriction to ofb{ti) ■ ■ ■ b{tm)v is a common 
eigenfunction of the operators T{w) £ End(Fun^(T^)[0]). 

5. Interaction-round-a-face models 

In this section, we consider a special case of the above construction, and 
relate T(w) to the transfer matrix of the interaction-round-a-face (IRF) (also 
called solid-on-solid) models of [Q, [ABF|. Therefore, our formulae give in 



particular eigenvectors for transfer matrices of IRF models. 

Notice first (see 0) that if we defin^a "Boltzmann weight" w{a, b, c, d; z), 
depending on complex parameters a, b, c, d, z, such that a — b,b — c, c — d,a — 
dG {1,-1}, by 

R(z, —2r]d)e[c — d] e[b — c] = w{a, b, c, d; z)e[b — a] e[a — d], 

a 

(the sum is over one or two allowed values of a) then the dynamical quantum 
Yang-Baxter equation translates into the star-triangle relation 

J2g w{a, b, g,f;z- w)w{f, g, d, e; z)w{g, b, c, d; w) 
= Eg w{f, a, g, e; w)w{a, b, c, g; z)w{g, c, d,e;z - w). 

We let W = be the tensor product of fundamental representations with 
evaluation points zi, . . . ,Zn- Then 

n 

L{z, A) = R^°^\z-zi, X-2r] ^ /i^^')) • • • i?(0'"-i)(z-z„_i, A-2?7/i("))i?(°")(z-z 

i=2 

In this formula the factors oiVmV(^W = y^i^+i are numbered from to 
n. The module is a highest weight module with highest weight functions 
of the form (Q) with pj = zj, qj = Zj + 27]. 

Let us now introduce a basis |ai, . . . ,a„) of Fun^(VK[0]), labeled by Oj G 
/X + Z with ai — at+i G {1, —1}, i = 1, . . . , n — 1, and — ai G {1, —1}- We 
let (5(A) = 1 if A = and otherwise. Then we define 

|ai, . . . ,a„)(A) = 5{X + 2ryai)e[ai - 02] e[a2 - a^] (g) ■ ■ ■ (g) e[an - ai]. 



^We adopt here a convention for the definition of the Boltzmann weights which is 
shghtly different than in |F[| and in better agreement with the literature. 
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Figure 1 . Graphical representation of the row-to-row trans- 
fer matrix of an IRF model. Each crossing represents a Boltz- 
mann weight w whose arguments are the labels of the adjoin- 
ing regions and the difference of the parameters associated 
to the lines. 



If r is the shift operator r/(A) = /(A — 27]), then T\ai, . . . , an) = \ai — 
l,...,a„ — 1). Using this, and the fact that h^^^ai, . . . ,0^) = (aj+i — 
Oj)|oi, . . . ,a„), we get 

n 

T{z)\ai,...,an) = ^ Y['^(.bj,aj,aj+i,bj+i;z-Zj)\bi,...,bn), 

bi,...,b„ j=l 

with the understanding that bn+i = bi, an+i = ai. The (finite) sum is over 
the values of the indices bi for which the Boltzmann weights are defined. 
Comparing with [Q, we see that T(z), in this basis, is the row-to-row transfer 
matrix of the (inhomogeneous) interaction-round-a-face model associated 
to the solution w{a,b,c,d; z) of the star-triangle equation (see The 
situation is best visualized by looking at the graphical representation of Fig. 

s- 

This construction can be in principle extended to higher representations, 
and we obtain in this way eigenvectors of transfer matrices of the IRF models 
of [pjMQj. 



6. The eight-vertex model 

We show in this section how to obtain from our result eigenvectors for 
the transfer matrix of the eight-vertex model. 

The eight-vertex model is based on Baxter's solution 

Rsviz) = asv{z){Ei^i ^1,1 + ^-1,-1 ® ^-1,-1) 
+bsviz)iEi^.i + 
+C8v{z){Ei^-i ^ ^-1,1 + ^-1,1 ^1,-1) 
+dsviz){Ei^i ^-1,-1 + ^-1,-1 ^ ^1,1), 
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of the Yang-Baxter equation. The coefficients are 

a8V{z) = — , bsv[z) — 



9i{z-27])9o{0)' '''' ' 9o{z-2r])9o{0)' 
in terms of the theta functions with characteristics 

Remarks. Baxter uses —z instead of z. The classical Jacobi notation, used 
by Baxter, is 0i{z) = H(2Kz), ^o(^) = e(2i^z). 

For any (generic) zi, . . . ,Zn, one then defines commuting transfer matrices 
Tsv{z) = tvoRsviz - -2i)(°^) ■ • • Rsviz - ^n)(°"^. (5) 

acting on (C^)®". The relation with our transfer matrices is based on the 
following identity, which is a version of Baxter's vertex-IRF transformation: 

Proposition 7. Let S{z, A) be the matrix 

_J_f9o{z-X + l/2) -9o{-z-X + l/2y 
9{X) \-9i{z - A + 1/2) 9i{-z - A + 1/2) 

Then 

S{w, A)(2)5(z, A-2r//i(2))(i)i?(2-u;, A) = Rsv{z-w)S{z, X)^^'^S{w, X-2r]h^^^)^'^^ 

Proof: We need to recall some well-known facts about the functions 9 a and 
Baxter's i?-matrix R^v- 

The functions 0^, a = 0, 1 arc entire, and uniquely determined up to nor- 
malization by the properties 0q(z+1) = {—lY9a{z), Oa{z+T) = ie~''*(^+'^/^)6'i_, 
Also, 9a{—z) = {—l)°'9a{z), and the zeros of 9a are simple and of the 
form r + 2st if a = 1 and r -|- (2s -|- l)r if a = 0, (r, s G Z). We have 
9o{z)9i{z) = C{t)9{z) for some constant C(r), as one can see comparing 
transformation properties under translations by Z + rZ. 

Let A = i }^ \ \, B = \ Then R^v{z) commutes with A® 



^0 ly vi 

^ and 5 B. Moreover Rw{z + 1) = A^^'> Rsv{z)A^^\ Rsv{z + r) = 
g-27ri7,^(i)^g^(^^^)^(i)^ and RsviO) = P, the Qip u (g) v ^ v (gi u. As a 
function of z, Rsv is meromorphic. Its poles are simple and are at 2ri 
modulo Z -I- rZ. The residue of Rgviz) at z = 2r] is 9{2r])/9'{0) times the 
anstisymmetrization operator Il:ufi^v^ufi^v — vfgu. 

Let (f){z) = {9i{z),9o{z)). This vector is the (up to normalization) unique 
vector with entire holomorphic components such that (l){z + 1) = A(j){z), 
^(z + r) = ie-^'^(^+^/2)^^(^)_ 

We claim that the statement of the Proposition is equivalent to the fol- 
lowing set of identities, which, incidentally, is essentially the standard form 
of the vertex-IRF transformation. 
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Lemma 8. Let (j) {z, A) = ipi'^z — A + ^). Then 

Rsviz-w)(l)^{z,X^27])^(j)^{w,X) = (/)=^(z,A) (/)^(t«,AT2r/) 

Rsviz - w)(l)^{z,X±2r]) ^ (j)'^{w,X) = a{z - w,±X)cl)^{z, X) c/)'^ {w, X^ 2r]) 

+f3{z - w, ±A)</>=F(z, A) (g) (j)^{w, X±2r]). 

The (known) proof of this lemma consists in comparing transformation 
properties under lattice translations and poles of both sides of the equations 
as functions of z, w, and using the uniqueness of (j). 

It remains to show that these identities are equivalent to the statement 
of the Proposition. Let us write them in matrix form: Let S{z, A) be the 2 
by 2 matrix whose columns are (j)'^ and then the previous lemma reads 

R8viz-w)S{w,X)^^^S{z,X-27jh^^^)^^^ = S(z,A)«^(u;,A-2??/iW)(2)i?( z-w,X) . 

The transposed of a matrix in End(C2oC2) is defined by the rule (X^y)* = 
X* y*, X,Y £ End(C^). The statement of the Proposition is proved 
by transposing both sides of the equation, and inverting the matrices S. 
The identity we need is = const S, which follows from the identity 

det S{z, A) = const 9{z)6{X). The latter formula can be proved by comparing 
the transformation properties under lattice translations of z and A, and 
using the fact that 9 is the unique function (up to normalization) such that 
e{z + 1) = -e{z) and e{z + r) = -e-^^(^+2^)0(z). □ 

The transfer matrix T{z) of the previous section has the form T(z)^(A) = 
a(z, X)i){X - 277) + d{z, X)i){X + 277), where a(z, A), d{z, X) £ End((C2)^"[0]) 
are the diagonal matrix elements of L{z, A). 

Corollary 9. Fix generic zi, . . . ,Zn G C. Let 

SniX) = S{z^, A)(")S(z„_i, A - 277/i(-))("-i) • • • 5(zi, A - 2,? h^^^fl 

Then, for all generic complex z, X, 

Tsviz)Sn{X) = Sn{X + 2r/)a(z, X + 2r]) + 5„(A - 2??)d(z, A - 2r/), 
on the zero weight space (C^)®"[0]. 

Proof: Let us write T^vi^) = tro-^^sy (^)) see (^). The matrix Lgy^z) acts on 
the tensor product of n + 1 copies of C^, numbered from to n. By iterating 
Proposition we obtain 

Sn{X)S{z, X - 27]hf ^L{z, A) = Lsv{z)S{z, A)W5„(A - 27?/i(o)). (6) 

In this formula Sn{X) acts on the factors numbered from 1 to n, and h = 
hW + . . . + /iW. Let L^{z, A) be defined by 

L{z,X)e[a] = e[/3] L^(z, A), 
/3=±1 
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and set t/ja = S(z, A + 2T/a)e[a]. Replacing A by A + 2rja in (^) and acting 
on a vector of the form e[a] ®u, a & {1, —1}, where hu = 0, yields 

^^p® Sn{X + 2r]a)L'^{z,X + 2r]a)u = L^v{z){'ipa O 5'„(A)n), 

where we used the fact that L{z, X) commutes with /i'-'^-' + h. Since, for 
generic parameters, form a basis of C^, and L\{z,X) = a(z,A), 

LZ\{z, A) = d{z, A), the proof is complete. □ 



Theorem 10. Let f ^ J /(A) be a linear function on a space T of functions 
of X G C, such that f f{X + 2r]) = j /(A) for all f G T. Extend j to vector- 
valued functions by acting componentwise. Then for each eigenfunction V'(A) 
ofT{z), the vector J S'„(A)^(A), if defined, is an eigenvector ofTsv{z) with 
the same eigenvalue. 

This theorem is an easy consequence of the previous corollary. What is 
left to do is to find linear forms J defined on the components of S{X)'ip{X) 
for the Bethe ansatz eigenfunctions ^ of Theorems ^. 

This can be done easily in two situations: notice that both ^(A) and ip{X) 
are periodic in A with period 2, if the parameter c belongs to vriZ. 1. If 
rj = p/q is rational, we may choose 

/ /(A) = j;/(/. + 2r?i). 



for any generic //. 
2. If is real, we set 



for generic /x. 



j=0 



7. g-DEFORMED LAME EQUATION 

We consider here the special case where W is the evaluation module 
V2m(0) with positive integer m. Then the zero weight space is one-dimensional. 
From the expression for a and d given in |FV3|, we see that the transfer ma- 



trix has the form (cf. Q , where a similar observation is made in the context 
of the Sklyanin algebra) 

^^'^ = e{z-%mti)v)^' 

where the difference operator L is independent of z and is given by 
^(A + 27?m),., ^ , , 9{X-27]m) ^ ^ ^ 

' " — e{x) — ~ ^' — e{x) — 
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This operator is in fact conjugated to the two-body Ruijsenaars operator 
[Q. It obeys {L(l),il)) = {(p,Ltlj) with respect to the symmetric bihnear form 

^<p,^p) J njL^0(A_2r?j)0(A + 2r?j)' 

where / is a linear function invariant under shift by 2rj, and change of sign 
of the argument, defined on a suitable space of functions. If rj is real, we 
may take / to be the integral over 7 = 7+ + 7_ , where the straight path 7+ 
goes from r/2 to 1 + t/2 and 7_ goes from —1 — t/2 to — t/2. This linear 
form is defined, say, on the space of meromorphic functions whose poles are 
contained in {2'qj + k + lT,j, k, I G Z}, which is preserved by the action of 
L. 

Note that L has periodic coefficients, and therefore preserves the space 
of Bloch functions ■0 such that + 1) = ^tp{\). The bilinear form is 
well-defined on this space of functions. 

We consider the eigenvalue problem 

Lil; = €-0. (8) 

It is a g-deformation of the generalized Lame equation: we have the expan- 
sion of L in powers of i]: 

The differential operator appearing in the second order coefficient is the 
generalized Lame operator (up to conjugation by ^(A)™). Our Bethe ansatz 
solution is a g-deformation of Hermite's solution of the (generalized) Lame 
equation (see the last pages of [ [WW] ]). 

Theorem 11. Let (ti, . . . ,tm,c) be a solution of the Bethe ansatz equations: 
9{t, - r?(l + 2m)) e{t, - - 2^) _ 
9{t,-7j{l-2m)) ^^^eit,-t, + 27^) ' 

such that ti 7^ tj mod Z + rZ if i ^ j. Then 

m 

H^) = e'^l[e{\ + tj-v), (10) 

is a solution of the q-deformed Lame equation Lip = etp with eigenvalue 

^ ^ ^_2^ J{Aiim) ^ d{t, + (2m - 3)r?) 

e{2r]m) e{tj + (2m - l)r/) ' 

Proof: This theorem is a special case of Theorem ^. It follows from (0) 
that the formula for the eigenvalue is e = ~^~^^^3^y^£('2), where e{z) is 
the function in Theorem ^. Since this expression is independent of z, we 
can evaluate it at any z. The formula given in this theorem is obtained by 
taking z = {1 — 2m)r], so that the second term in e{z) vanishes. □ 
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The Bethe ansatz equations have the form bi{t) = e'^^'^, i = 1, . . . , m. We 
may ehminate c and consider the set of {ti, . . . ,tm) obeying the equations 
bi{t)bj{t)~^ = 1. The functions are doubly periodic meromorphic 

functions with periods 1 and r in each of the variables tj. Also, eigenfunc- 
tions corresponding to solutions related by shifts of the variables ti by 1 or 
r are proportional. Therefore the set 

X = {te (C/Z + TZr/S^\bi{t) = bj{t) ^oo,i = l,...,m, tj{i + j)} 

is an algebraic subvariety of the symmetric power of our elliptic curve, which 
parametrizes our eigenf unctions. We call it the Hermite-Bethe variety. On 
this variety, we have a single- valued function e^, the square of the eigenvalue, 
which is the restriction to X of a doubly periodic function. Therefore it is 
an algebraic function : X — > C. The eigenvalue, which is the square 
root of this function, is two- valued on X. This is connected to the fact that 
generically there are two eigenfunctions corresponding to a single point of the 
Hermite-Bethe variety: a point of X determines c only modulo (4r/)~^27riZ. 
If is an eigenfunction corresponding to a point of X, then e'^^^^'^^ilj{X) 
is an eigenfunction with the opposite eigenvalue corresponding to the same 
point but with c translated by 27ri/4ry. 

Note that the space of meromorphic solutions of the difference equation 
(^) is a vector space over the field K of meromorphic 2r/-periodic functions 
of A. If (ti,...,t m, c) is a solution of the Bethe ansatz equations, then, for 
all /c G Z, {ti, . . . ,tm,c + 2iTik/2rj) is also a solution with the same value 
of e. The corresponding eigenfunction is proportional (with a 2?7-periodic 
coefficient) to the original eigenfunction. Therefore we should consider c 
modulo nirj^^'Z. 

Note also that we have an involutive automorphism a of X, sending 
(ii, . . . , tm) to (2?7 — ti, . . . ,2r] — tm)- If {t, c) is a solution of the Bethe 
ansatz equations, then also (cj(t),— c). The corresponding eigenfunctions 
are related by the "Weyl reflection" ^ 'ip{—X) and have the same 

eigenvalue. 

We now turn to the question of completeness. Let S'^Ej- = {'L + t'L)'^ / Sm 
be the symmetric power of the elliptic curve. 

Theorem 12. Suppose that rj ^ C is generic. For generic e G C, there are 
precisely two solutions 

{ti, . . . ,tm,c) and {2r] - ti, . . . ,2r] - tm,-c), 

in S"^Et- X (C/-n'iri~^7j) of the Bethe ansatz equations (|^) with given e. The 
corresponding eigenfunctions ip± are linearly independent over the field K 
of 2r]-periodic meromorphic functions of \, and all solutions of the q-Lame 
equation ^ are linear combinations of tpj^, "0-- with coefficients in K. 

Proof: It is straightforward to check that the reflection ti i— > 2ry — tj, c i— > — c, 
maps solutions to solutions preserving the value of the eigenvalue e. 

Let X be the closure of X in the symmetric power of the elliptic curve. 
Then X contains the point P = {{l — 2m)r}, . . . , —2>r}, —rj). In a neighborhood 
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of P we can introduce local coordinates ui = ti — (1 — 2m)rj, uj = tj — tj-i — 
'^Tjj 2 < J < m. In these coordinates, X is described by the equations uj = 
CjUm + 0{um)), i = 1, . . . , m — 1, for some constants cj. Therefore, in 

this neighborhood, X is a non-singular curve, and local parameter 

at P. The eigenvalue has the form 

e = const ^"^/^(l + 0{um)) 

in a neighborhood of P. Therefore e is a non-constant function on X. Sim- 
ilarly one shows that c appearing in the Bethe ansatz equation is a non- 
constant function on X (c diverges at P). Since is algebraic, we have, for 
any generic value of e, a solution (ti, . . . , tm, c) of the Bethe ansatz equations, 
and thus an eigenfunction tp^. Let V'- be the eigenfunction with the same 
eigenvalue associated to the reflected solution (2ry— ti, . . . , 2r]—tm, — c). Since 
generically ip^ and ip- have different multipliers, they are linearly indepen- 
dent over K: suppose ■0+ = Cip- for some C € K. Then C is 2?7-periodic 
and obeys C(A-|- 1) = e^'^C(A). If rj is real and irrational and e^'^ 7^ 1, we get 
a contradiction, so '0+ and -0- are linearly independent. The case of generic 
r] is treated by analytic continuation. 

Next we show that every solution is a linear combination of tp^ and tp-. 
The (difference) Wronskian W{f,g) of two functions /, g is the function 
/(A -|- 2r])g{X) — f{X)g{X + 2r]). Two meromorphic functions are linearly 
dependent over K if and only if their Wronskian vanishes. If /, g are solutions 
of (P) then their Wronskian obeys the difference equation W{X + 2r]) = 
u{X)W{X), where the function u is a combination of the coefficients of the (^. 
Thus for any solution /, the functions A± = ±W{f, ipz^)/W{'ip+,tp-) are 2r]- 
periodic. On the other hand, A±{X) are (by Cramer's rule) the coefficients 
in the expression of (/(A), /(A -|- 2r])) as a linear combination of the linearly 
independent vectors {'ip±iX),'ip±{X + 2r])). In particular, 

/(A) = A+(A)V^+(A)+^_(A)V^_(A). 

Therefore / is a linear combination of ip± with 2r7-periodic coefficients. 

It remains to prove that for each generic e there are not more than two 
solutions of the Bethe ansatz equation with given e. Suppose that there were 
a third solution {t[, . . . ,t'^,c') distinct from the two we have constructed. 
In particular {t'l, . . . ,t'^) represents a point in the symmetric power of the 
elliptic curve which is distinct from the two points represented by (ti, • • • > ^m) 
and {2r] — ti, . . . ,2r] — tm)- Let ^p' be the corresponding eigenfunction. Thus 
tp' = atp-\- + bip- with 2?7-periodic coefficients o, b. We have that a = 
W(^p' ,'^p-)/W{^p+,^p-). The three functions ip' ^tpij^^ip- have all the form 
(|To|). Therefore the meromorphic function a obeys the equations o(A + 2r/) = 
o(A), o(A + 1) = Cia(A), a(A + t) = C2a{X) for some constants Ci,C2- If 
rj is generic, this implies that a is of the form a^e^^^/^ for some constant 
oo and some integer s (a cannot have poles since it would have a dense set 
of singularities). Similarly, b has the same form. We may moreover assume 
that e'^ is generic, so comparing multipliers we see that either a = or 6 = 0. 
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But the zeros of ip' are not equal to the zeros of "0+ or to the zeros of ^|J- , a 
contradiction. □ 

We see from the proof of this result that the irreducible component(s) of 
X containing P and its reflected point cr(P) form a curve Y C X. On this 
part of X, takes every generic value precisely twice. 

Let us summarize our results. 

Theorem 13. The closure X in S'^Er of the Hermite-Bethe variety con- 
tains an algebraic curve Y. It is a two- fold ramified covering : y ^ of 
the Riemann sphere. It has an involutive automorphism a : Y ^ Y permut- 
ing the sheets and preserving e^. For each generic point t £ Y, there are two 
solutions {t,c), {t,c-\-iTr/2r]) of the Bethe ansatz equation inYxC/TTir]~^7j. 
The corresponding eigenf unctions are related by ip{X) e'^'^'^/'^'>l>{\) , and the 
eigenvalues are the two square roots of e'^{t). 

One way to formulate this result is that eigenfunctions are parametrized 
by the spectral curve, the double covering of Y on which e is single-valued. 

8. Proofs 

This section contains the proofs of Lemma ^ and Theorem ||. These proofs 
are based on the following technical result. 

Lemma 14. Let V{z) with basis e[l],e[— 1] denote the fundamental repre- 
sentation of E^-^riish) with evaluation point z (see Sect. |^. Let ti, . . . ,tm, 
zi, . . . , Zm, T,ri be generic complex numbers, such that Im(T) > 0, and let 
W = V{zi) • • • V{zm). Then the 2"" vectors 

([lKtj))i4M^---^4M) e Fun(T^), 

where J runs over all subsets o/ {!,..., m}, are linearly independent over 
the field of meromorphic functions of A. 

Proof: Let us denote these vector by wj. It is sufficient to show that the 
values wj{\) are linearly independent for some value of A, and some value 
of the parameters. This will be shown by considering the matrix relating 
wj{X) to the basis ej = e[ai] (8) • • • e[am] with aj = 1 iff j G J, and 
showing that, in some limit of the parameters, this matrix is upper triangular 
with respect to the lexicographical ordering of binary numbers: J > K iS 
> The limit is obtained by first taking r ioo which 

amounts to replacing 9{x) by sin(7rx), then A — > oo. In this limit, b{t) acts 
on V{z) as b{t)e[-l] = 0, b{t)e[l] = B{rj){l - exp(-27rf (z - 1 + 2rj)))-^ e[-l] 
for some constant B(r]) ^ 0. Now let us choose tj = Zj + e for some fixed 
generic e, and set Zj = \/—lNj, j = 1, . . . ,m. Let e[a]k be basis vectors 
of V{zk). Then, as — > oo, the vector b{tj)e[l]k tends to zero if j < A: 
and tends to a nozero multiple of e[— 1]^ if j = k. On the other hand, the 



ALGEBRAIC BETHE ANSATZ 



17 



matrix elements of a and d tend to non-zero finite values in this limit. But 
the action of b{t) on a tensor product ui • • • (8) Um consists of m terms, the 
j'th one being given by the action of b on the jth factor and the action of a 
or d on the other m — 1 factors (see |FV3|). 

It follows that • • • b{tj^) is a linear combination of vectors ex where 
K = {ki, . . . ,kr} with ki > ji, ■ ■ ■ ,kr > jr, and that the diagonal matrix 
elements (such that K = {ji, . . . , jV}) are non-zero. □ 



Proof of Lemma |^. The coefficients Aj, Dj can be in principle computed 
using the commutation relations repeatedly, giving them as some universal 
polynomials in the values of the functions a, (3, independent of the high- 
est weight modules we are considering. Since the 6's commute, we ob- 
tain m\ ways of representing a{w)b{ti) ■ ■ ■ b{tm) as a linear combination of 
b{ti) ■ ■ ■ b{tj^i)b{w)b{tj+i) ■ ■ ■ b{tm)a{tj) (and similarly for d). Since by the 
previous Lemma there exists a vector in a module so that these operators 
applied to this vector yield linearly independent vectors, it follows that all 
these representations must coincide. □ 



Proof of Theorem |^. We have to write an explicit formula for 6(ti) . . . b{tm)vQ, 
where vq is a tensor product of highest weight vectors of some highest weight 
modules. Let us first consider the case where vq = n^w is the tensor product 
of two highest weight vectors. Then, using the rules for the action on ten- 
sor products (see | FV3|| ) and the commutation relations, we obtain a linear 



combination of terms of the form 

T{-2r]h^^^)[b{U,)- ■ ■ biUMtj,) ■ ■ ■ a{tj^_J]umtk,)-bitkr^-Mith)-ditJv, 

where I = {ii < ... < ij, J = {ji < ... < jms}, lU J = {1,2, ...,m}, 
K = {ki < ... < A;™_,}, L = {h < ... < Is}, KUL = {l,2,...,m}. 
The notation T{—2r]h^'^^) indicates that the argument of the matrix valued 
function of A in the square bracket must be shifted by —2r] times the weight 
of the vector in the second factor. 

The coefficients can be computed using only the commutation relations 
and are therefore independent of the choice of highest weight modules. It is 
convenient to take u and v to be vectors of the form of the previous Lemma, 
since it then follows from the linear independence that the coefficients are 
uniquely determined. 

Suppose that we compute the coefficients by first applying b{tm) to the 
tensor product, then b{tm-i) and so on, and then use the commutation 
relations to shift d's and a's to the right of the Vs. Then it is clear that the 
coefficients of the terms with I B 1 and K B 1 vanishes. Similarly, if we use 
the fact that the 6's commute to act with b{tj) at the end, we see that the 
coefficients of the terms with I D K 3 j must vanish. We conclude that the 
only terms appearing with non- vanishing coefficients must have I D K = 9. 
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Therefore 



llb{U)l[a{tj] 



u0Ylb{tj)Y[d{ti 

30 «e/ 



for some universal, uniquely defined, coefficients Cj(ti, . . . , tm), / C {1, . . . , m}. 
By the commutativity of the 6's, these coefficients can be computed in dif- 
ferent ways with the result that 

for all permutations a £ Sm- Therefore, it is sufficient to calculate C/ for / = 
{1, . . . , s}, for all s E {1, . . . , m}, which can be done straightforwardly using 
the tensor product and commutation rules. Taking u equal to the highest 
weight vector of an evaluation representation Va(-z) and v a tensor product 
of such highest weight vectors gives a recursive procedure to compute all 
coefficients given in Theorem ^ □ 
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